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Recent experiments of the quasi-one-dimensional spin-i antiferromagnet Copper Benzoate es- 
tablished the existence of a magnetic field induced gap. The observed neutron scattering intensity 
exhibits resolution limited peaks at both the antiferromagnetic wave number and at incommensurate 
wave numbers related to the applied magnetic field. We determine the ratio of spectral weights of 
these peaks within the framework of a low-energy effective field theory description of the problem. 
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I. INTRODUCTION 



Recent experimentsi*^ have investigated the behaviour 
of the quasi-one-dimensional spin-1/2 antiferromagnet 
Copper Benzoate in a magnetic field. Neutron scattering 
experiments^ established the existence of field-dependent 
incommensurate low-energy modes in addition to low- 
energy modes at the antiferromagnetic wave number. 
The incommensurability was found to be consistent with 
the one predicted by the exact solution of the Heisenberg 
model in a magnetic field. However, the system exhibited 
an unexpected excitation gap, induced by the applied 
field. A theory for this effect was put forward by Os- 
hikawa and Affleck in Ref . Q . Application of a uniform 
magnetic field H induces a staggered field perpendicular 
to H. The staggered field is generated both by a stag- 
gered g-tensor"* and a Dzyaloshinskii-Moriya (DM) inter- 
action. The same physical mechanism has been found 
also in other materials.4 The effective Hamiltonian de- 
scribing such field-induced gap systems is given by^ 



where 



h = -fH 



(1) 



(2) 



The constant 7 is given in terms of the staggered g- 
tensor^ and the DM interaction. For Copper Benzoate 
the exchange constant is J ~ 1.57 meV and the induced 
staggered field is much smaller than the applied uniform 
field, h<$:H. 

As long as /i <C J, or equivalently as long as the field 
induced gap A is much smaller than J, it is possible 
to describe the low-energy degrees of freedom of hamil 
in terms of a massive, relativistic quantum field theory. 
This low-energy effective theory is obtained by abelian 
bosonization and is given by a Sine-Gordon model with 
Hamiltonian densitji^ 
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[(a,$)2 + (d^ef] - ii{h) cos(/3e) 



(3) 



Here $ is a canonical Bose field, 9 is the dual field and 
the coupling (3 depends on the value of the apphed uni- 
form field and has been calculated in Refs. ^^6 7^8] by using 
the results of Ref. . The spin velocity v also depends 
on H and is shown in Fig. 9 of Ref. It is useful to 
define 



87r-^2 



(4) 



The spectrum of the Sine-Gordon model SGM in the 
relevant range of (3 consists of a soliton-antisoliton dou- 
blet and several soliton-antisoliton bound states called 
"breathers. "-Sill The soliton gap as a function of H and h 
was determined in Ref. @ in the regime A <g; -ff , where 
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with B = 0.422169. Equation deltaRG is applicable as 
long as H is sufficiently smaller than J or to be more pre- 
cise as long as the magnetization is small. For magnetic 
fields comparable to J it is better to use the following 
expressioE>iS 
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Here v — v/{Jao) is the "dimcnsionless spin velocity", 
ao is the lattice constant and c{H) is given below. The 
breather gaps are given by 



A„ = 2Asin ( -|- 



(7) 



II. DYNAMICAL STRUCTURE FACTOR 

The staggered/oscillating components of the spin op- 
erators are expressed in terms of the continuum fields 
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FIG. 1: Schematic structure of the lowest-energy excited 
states relevant to neutron scattering experiments. Soliton 
and antisoliton occur in the vicinity of the incommensurate 
wave numbers tt ± 25 and are seen in S^^ , whereas the first 
breather Bi occurs in the vicinity of tt and contributes to 
S^^ . At higher energies further breather bound states as well 
as multi-particle scattering continua are present. 



and 9 as 

^ (-1)" c{H) cos{(3 6) , 
Sy ^ i-ir c{H) sin(/3 6) . (8) 

Here x = nao and the incommcnsuration S is determined 
from the exact solution of the Heisenberg model in a uni- 
form magnetic field^ [that is the Hamiltonian hamil for 
/i = 0]. The amplitudes a{H) and c{H) are at present 
not known analytically, but can be determined numeri- 
cally with high accuracy. We note that these amplitudes 
are also calculated in the absence of a staggered field, 
the expectation being that the changes due to a small 
h <^ H will be negligible. The data used in this work 
are obtained in the scheme of Refs. [TllIT^ : We calculate 
the spin polarization (S^) and the two-spin correlation 
function (S'^^S'J^/) in the Heisenberg chain of 200 spins 
using the density-matrix renormalization group method, 
and then, fit them to analytic formulas obtained from the 
abelian bosonization taking a{H) and c{H) as fitting pa- 
rameters. The results as well as other parameters, which 
are determined exactly^Sii are listed in Table for several 
typical values of the magnetization m. 

The inelastic neutron scattering intensity is propor- 
tional to 

/(c.,fe)«5:(l-^) (9) 
where a,P = x,y,z and the dynamical structure factor 



TABLE L Amplitudes a and c, the spin velocity v, the cou- 
pling /3, and the field H as functions of the magnetization 
m. The amplitudes are determined numerically except for 
m = 0.5 where exact values are shown. The figures in paren- 
thesis for a and c indicate the error on the last quoted digits. 



m 


a 


c 


V 




H 


0.02 


0.591(3) 


0.4937(3) 


1.54271 


2.35016 


0.17599 


0.04 


0.550(5) 


0.4883(2) 


1.51707 


2.31088 


0.34214 


0.06 


0.520(4) 


0.4863(2) 


1.48415 


2.27738 


0.50013 


0.08 


0.4947(6) 


0.4853(2) 


1.44425 


2.24653 


0.65001 


0.10 


0.475(1) 


0.4847(2) 


1.39796 


2.21731 


0.79164 


0.12 


0.454(1) 


0.4842(2) 


1.34593 


2.18927 


0.92489 


0.14 


0.437(2) 


0.4835(2) 


1.28879 


2.16216 


1.04965 


0.16 


0.422(2) 


0.4825(2) 


1.22720 


2.13587 


1.16589 


0.18 


0.4070(7) 


0.4810(2) 


1.16178 


2.11029 


1.27360 


0.20 


0.3938(8) 


0.4790(2) 


1.09314 


2.08538 


1.37287 


0.22 


0.3813(6) 


0.4764(2) 


1.02184 


2.06107 


1.46380 


0.24 


0.3700(8) 


0.4731(2) 


0.94844 


2.03735 


1.54656 


0.26 


0.3596(7) 


0.4690(2) 


0.87347 


2.01418 


1.62134 


0.28 


0.3499(4) 


0.4639(2) 


0.79741 


1.99153 


1.68839 


0.30 


0.3406(4) 


0.4578(2) 


0.72074 


1.96940 


1.74794 


0.32 


0.3330(2) 


0.4504(2) 


0.64387 


1.94775 


1.80030 


0.34 


0.3262(2) 


0.4416(2) 


0.56722 


1.92658 


1.84575 


0.36 


0.3200(3) 


0.4310(2) 


0.49116 


1.90586 


1.88462 


0.38 


0.3145(4) 


0.4183(2) 


0.41602 


1.88559 


1.91723 


0.40 


0.3094(2) 


0.4029(1) 


0.34212 


1.86574 


1.94390 


0.42 


0.3070(8) 


0.3841(1) 


0.26973 


1.84631 


1.96497 


0.44 


0.3058(2) 


0.3601(1) 


0.19912 


1.82727 


1.98079 


0.46 


0.3062(6) 


0.3284(1) 


0.13049 


1.80863 


1.99168 


0.48 


0.309(1) 


0.2802(1) 


0.06407 


1.79036 


1.99797 


0.50 


0.3183 








1.77245 
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S"'^ is defined by 

S^^{u;,k) = / -e-^'='-+-*(5r+i(t)5f(0)). 
1=1 •'-^ ^'^ 

(10) 

Here fc denotes the component of fc along the chain direc- 
tion. A schematic representation of which excited states 
will contribute to the various components of the dynam- 
ical structure factor at fc = tt/oo and = (tt ± 25) /a^ 
is shown in Fig. ^ At the antiferromagnetic wave num- 
ber there are several breather excitations and at higher 
energies multiparticle continua.^ These contribute to the 
XX and yy components of the dynamical structure fac- 
tor, which have been determined in detail in Ref. 0. 
At the incommensurate wave numbers A; = (tt ± 25) /oq 
there are soliton and antisoliton states and at higher en- 
ergies multiparticle scattering continua. In this paper we 
calculate the single particle soliton/ antisoliton contribu- 
tions to S^^ and compare them to the dominant feature 
in the dynamical structure factor, the contribution of the 
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lightest breather bound state Bi to S^^. where 
The hghtest breather Bi has a gap Ai and contributes 

to syy asi 



/ vr 
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Cy{H) = 2iJ c^{H) 



2cos«/2)\/2sin(7r^/2) exp 




' A V^ r((i + o/2) \ 
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X exp 



i V2sinh(/32t)sinh(87ri)cosh[(87r-/32)t] 47r 
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Here we have used the normalizations of Ref. (13| . The leading contributions to the longitudinal structure factor at 
the incommensurate wave numbers = (tt ± 25) j oq are due to soliton and antisoliton. Using the results of Ref. 0| 
we obtain 



where 



i5J 



X exp 



dt 



exp[-(l + Ot]-l 



2-K e 



-2t \ 1 



i V 2 sinh(^i) sinh[(l + ^)t] cosh(i) 2 sinh(t^) 
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Here the constants Ci,2 are given by 



Ci = exp ( — 



exp 4 



' dt smh^{t/2) sinh[t(^- 1)] 
Tsinh(2i) sinh(Ci) cosh(t) 
dtsmh^{t/2) smh[t{( - 1)] 



(15) 



t sinh(2t) sinh(^t) 

As was pointed out in Ref. 6], at _ff = the low-energy 
effective theory of hamil is SU(2) symmetric. In our no- 
tations this implies that 

(16) 



lim f",^ , 



— 2 lim „ , , . 



Equation su2 is easily verified numerically. In order to 
evaluate Cy^z we need to know the constant of propor- 
tionality 7 that relates the induced staggered field h with 
the applied uniform field H . This constant differs from 
compound to compound. On the other hand, 7 enters the 
expressions for Cy^z only via the soliton gap A. Hence 
it is useful to isolate the 7 dependence and consider the 
quantities 



A 



C'y{H) = Cy{H) - 



CiiH) = Cz{H) ( - 



(17) 



The amplitudes C'y ^ (H) are shown as functions of the 
magnetization in Fig. |21 



III. COPPER BENZOATE 

We are now in a position to determine the ratio be- 
tween the spectral weights of the first breather (seen in 
the transverse structure factor 3^^) and the soliton (seen 
in the longitudinal structure factor 5*^^). In an ideal situ- 
ation one would carry out measurements with momentum 
transfers only along the y and z directions respectively. 
In practice, the experiments on Copper Benzoate were 
carried out with momentum transfers fei.2 respectively, 
where 

(fci • a,fci • 6,fci • c) = 27r (-0.3,0,1) , 

(fca • a, fea • 6, fca • c) = 27r (-0.3, 0, 1.12) , (18) 

Here c points along the chain direction and the antiferro- 
magnetic wave number corresponds to 27r/c (c = 6.30^) 
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2: Amplitudes Cy^^{H) as functions of the magnetiza- 



as there are two copper atoms per unit cell along the c- 
axis. To make contact with our previous notations we 
need to set oq ~ c/2. The measurements with momen- 
tum transfers ki and fc2 probe the dynamical structure 
factor around tt/co and the incommensurate wave num- 
ber (tt -I- 25) /ao respectively. In order to make direct 
comparisons with the experiments we need to relate the 
(a, b, c) coordinate system describing the crystal axes to 
the {x,y,z) spin coordinates. By definition z and x are 
the directions of the uniform and staggered fields respec- 
tively. In the experiments of Ref. |3| the uniform field 
was applied along the 6-direction. Based on a polariza- 
tion analysis it was suggested in Ref. |^ that the stag- 
gered field lies in the ac plane and encloses an angle of 
a = —72° with the a-axis. In the vicinity of tt/oo the 
dominant contribution to the structure factor comes from 
the transverse correlators. This implies that 



J(w,fci) (X (0.083 cos^ a + 0.917 sin^a)^^^ ^,— 

V oo 

-|-(0.083 sin^ a + 0.917 cos^ a)S'''' ^w, ^ 

« 0.84 0.16 ( oj, —] .(19) 

On the other hand at momentum transfer k2 the dom- 
inant contribution to the structure factor is due to the 
longitudinal component 



/(C^,fe2) CX S^' 



25 



ao 



(20) 



As was pointed out in Ref. there are unresolved issues 
concerning the polarization analysis and the estimate of 
the angle a should be regarded with some caution. It is 
possible to infer a by analyzing other experiments such 
as specific heat and ESR measurements. The analysis of 
the specific heat data suggests that a « —82°^ which 
leads to a contribution of about 90% of S^^ in I. 



The neutron scattering experiments of Ref. |3| were 
performed in a uniform magnetic field of 7 T, which cor- 
responds to a magnetization per site of m « 0.06. The 
breather and soliton gaps were observed at 



A w 0.22 meV , Ai « 0.17 meV. 



(21) 



Using the expression deltaRG for the soliton gap we can 
infer the coefficient of proportionality between the uni- 
form field H and the staggered field /i as 7 « 0.06. 
Taking 7 — 0.06 and m = 0.06, Eq. deltaRG gives 
A = 0.215 meV, Ai = 0.171 meV and we will use this 
set of parameters for our further analysis. 

Under the above assumptions we may now determine 
the spectral weights of the coherent soliton and breather 
peaks in the dynamical structure factor. The results are 
shown in Fig. O 
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FIG. 3: Amplitudes Cy^z{H) as functions of the magnetiza- 
tion. 

At a magnetization of m = 0.06 we have 



a(i/) 



(22) 



The spectral weights are obtained by integrating the 
respective structure factors over frequency at fixed mo- 
mentum, i.e., 



duj syy uj 



TT 

TT + 2S 



(23) 



The ratio of spectral weights between the first breather 
Ibi and the soliton Ig is approximately 



Ib, _ Cy{H) A 



3.64. 



(24) 
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In order to compare to experiment, we need to take into 
account the different momentum transfers in the mea- 
surements of S*^^ and S'^^ respectively. Prom I and 12 we 
arrive at the following theoretical prediction for the ratio 
of intensities 

R = 0.84^ « 3.06. (25) 

The experimentally observec^ ratio of peak heights be- 
tween the breather and soliton peaks is approximately 
2.8. This is in reasonable agreement with our result. For 
a better comparison one should take into account the 
resolution function of the instrument in both momentum 



and energy, but this goes beyond the scope of our present 
analysis. 
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